FURTHER EXAMPLES WITH MOMENTS OF GAMMA 

TYPE 
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This is an appendix to 20(] containing further examples. See [20(] for 
notation and for examples and equations referred to below by numbers. See 
also the further references in [gO, Addendum]. 

This appendix will probably be extended with more examples in the fu- 
ture. 

Appendix B. Further examples 

Example B.l (Rayleigh distribution). The Rayleigh distribution R is 

the chi distribution with density xe~^ This is a special case of 

Example 3.6, and we have 

jg^s ^ 2*/2r(s/2 + 1), -2<Res<oo. (B.l) 

We have p+ = oo, p_ = -2, 7 = 7' = 1/2, 5 = 1/2, x = 0, Ci = vr^/^ 

Example B.2 (Maxwell distribution). The Maxwell distribution M is 
the chi distribution x(3), with density (2/7r)^/^x^e~^ This is a another 
special case of Example 3.6, and we have 

ns/2 . os/2+1 , o 

^^' = fk'(^^)=^^(l + 5)- -3<Re.<0.. (B.2) 

We have p+ = 00, p_ = -3, 7 = 7' = 1/2, 5 = 1, x = 0, Ci = \/2. 

Example B.3 (Type-2 Beta distribution). The type-2 Beta distribution 
[341 . Chapter 4] has density 

/(x) = ^^" + ^^ x"-Hl + x)-("+^\ x>0, (B.3) 

r(a)r(/3) 

for two parameters a, /3 > 0. A variable X^^p with this distribution has 
moments given by 

A comparison with (3.1) shows that Xa^p = Ta/T'^, with and F^ inde- 
pendent. In particular, see Example 3.7, the F distribution is of this type 
(up to a constant factor): Fm,n = (^^/^)-^m/2,n/2- 
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nXo.,p + 1)^ = Ei3^7^ = ^)JJ> Res</3. (B.5) 



We have /0+ = /3, /9_ = —a, 7 = 2, 7' = 0, 5 = a + /3 — 1, x = 0, 
Ci = 27r/(r(a)r(/3)). 

Note that 1 + = S^^^, the inverse of a (usual) Beta distributed 

variable, see Example 3.4. Thus l + Xa^p also has moments of Gamma type, 
with 

r(a + /3)r(/3 
r(/3)r(a + /3 

Example B.4 (Cauchy distribution). The Cauchy distribution with den- 
sity l/(7r(l + x^)), —00 < X < 00, equals the t-distribution in Example 3.8 
with n = 1. Hence, if X is a random variable with a Cauchy distribution, 

then |X| = 1 71 1 = F^^^ and |X| has moments of Gamma type 

E|X|" = -rf- + -Vf---') = 7^-^, -l<Res<l. (B.6) 

' ' vr V2 2/ V2 2/ cos(7rs/2)' ^ ' 

Cf. Example 3.19, where A = — log \X\. 

We have p+ = 1, p_ = -1, 7 = 1, 7' = 0, 5 = 0, x = 0, Ci = 2. 

Example B.5 (Beta product distribution). Dufresne [ll[ has shown 
that if a, 6, c, d are real, then there exists a probability distribution G(a, c; a+ 
6, c + c?) on (0, 1) with moments 

vx^- r r(a + g)r(c + ^) 

-^r(a + 6 + .)r(c + d + s) ^^-^^ 

(where necessarily C = T{a + &)r(c + d) / (T{a)T{c)))^ if and only if either 

(i) a > 0, c > 0, 6 + (i > and min(a + 6, c + d) > min(a, c), or 

(ii) (IBTt]) degenerates to EX* = CT{a + s)/r(a + /3 + s) with q > 
and /3 > 0, so X has a Beta distribution or X = 1. (This degenerate 
case occurs if 6 = 0, d = 0, a + 6 = c or c + d = a.) 

We have /)+ = 00, p- = — min{a, c}, 7 = 7' = 0, 5 = — 6 — d, x = 0, 
Ci = C. 

The case when all a,b,c,d > is just a product of two independent 
Beta variables Ba,bBc,d, see Example 3.4, but there are also other possible 
parameter values, for example (2,5,8,-1) given in [ll| . 

Moreover, if we allow complex parameters a, 6, c, d, there is exactly one 
more case [ul], viz. a>0, c>0, 6 + d>0 (entailing Im(6) = — Im(d)), and 
Re(a + h) = Re(c + d). In particular, we may take a = c > and d = b 
for any complex b with Re 6 > 0. However, complex parameters are not 
included in the class of distributions studied in see Remark 11.3. 

Example B.6 (Density of ISE). The ISE (integrated superbrownian ex- 
cursion) is a random probability measure introduced by Aldous It was 
shown in [6] that the ISE a.s. is absolutely continuous, and thus has a (ran- 
dom) density /ise(x), x £ (—00,00). 

The ISE can be described as the occupation measure of the head of the 
Brownian snake, see Le Gall [H, Chapter IV] or Le Gall and Weill [sO] for 
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details; see also [19|, Section 4.1]. Thus /ise(^c) is the local time of the head 
of the Brownian snake. Moreover, fis^ix) arises for example as a limit of the 
vertical profile of random trees, see [33], 0], 0], [1] and [lol ]. 

The distribution of /isE(a^) for a fixed x is given by a rather compli- 
cated formula, see 0] and 0|; in the case x = it simplifies and /ise(O) = 
21/42-1^^^1/2^ where ^2/3 is a positive 2/3-stable variable with Laplace trans- 
form E e-^-^^/a = g-*"/ 

, see Example 3.10. Thus /ise(O) has moments of 

Gamma type with 

EU^iOy = 2^/43-^ r(3s/4 + l) < Res < 00, (B.8) 

r(s/2 -I- 1) 

see [i] and 0]. We have p+ = 00, p- = -4/3, 7 = 7' = 1/4, 5 = 0, 
x=-flog2-ilog3, Ci = y372. 

Example B.7 (Average ISE). The ISE in Example lB.Gl is a random prob- 
ability measure ^ise! taking the expectation we obtain a deterministic prob- 
ability measure E /Xise; which is the distribution of a random variable X that 
can be seen as a random point given by a random ISE. (This is, for example, 
the limit distribution of the label of a random node in a random tree under 
suitable assumptions and normalizations.) X has a symmetric distribution, 
and \X\ has moments of Gamma type with 



23^/4 /s 1 



s 



ElXl^ = ^r^- + -jr^- + l), -l<Res<oo, (B.9) 



see 0] and 0. We have p+ = 00, p_ = -1, 7 = 7' = 3/4, S = 1/2, 
x= -|log2, Ci = ^. 

Example B.8 (Blocks in a Stirling permutation). Let A; > 2 be a fixed 
integer. It is shown in ^] that the number of blocks in a random ^-Stirling 
permutation of order n (see [12] for definitions) after suitable normalization 
converges in distribution as n — )• cxd to a random variable C, with moments 
of Gamma type given by 

r(i+|) _^(^ , 1^ r(s + 2) 

k 



EC' = (s + l)!— ^ ff- = r 1 + - / -2<Res<oo. 



(B.IO) 

As explained in |22| . this is actually a special case of (9.1). 

We have /3+ = 00, p- = —2, 7 = 7' = (5 = (A; — 1)/A;, x = ^logA:, 
C7i = A:('=+2)/2fer((A: + l)/A;). 

Example B.9 (Distances in a sphere). Let X\ and X2 be two indepen- 
dent random points, uniformly distributed in an n-dimensional ball of radius 
a, and let D := \X\ — A2I be the distance between them. (Here n > 1.) 
Note that < D < 2a, so D/2a G [0, 1]. Hammersley [ll] showed that the 
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density function oi D/2a is 



/„(A) = 5%±il^ /V d= (B,U) 



A 



and as a consequence, for Res > — n, 

^(2^+2) (n + s)r(n+ + 1) 

and, equivalently, 

r(s + n)r(is + in + i) 



ED" = C7(2a)' 



r(s + n+l)r(is+n+l) 



C'a'^, 1 ^^^+'^1 (B.13) 



with C = nT{n + l)/T{\n + \) and C" = 7ri/22-"C. (Hammersley \il 
did not specify the range of s, and presumably intended only positive and 
perhaps integer values, but the formula follows by (jB.lip for any s with 
Res > n. Alternatively, the result extends from positive s by Theorem 2.1.) 

D thus has moments of Gamma type, with = +00, p_ = — n, 7 = 
y = 0, ,5 = -{n + 3)/2, X = log(2a), Ci = 2("+i)/2C. 

It follows from (1B.13|1 that if Bn^i and i?(n+i)/2,(n+i)/2 are independent 

Beta distributed variables, then E(D/2a)'^ = Ei?^ i-^(*n+i)/2 (n+i)/2' 
ample 3.4, and thus 

^ = 2«5n,i<+l)/2,(n+l)/2- (^-14) 

Cf. Remark 1.5. Note further that Bn,i = C/^/" where U ~ U(0, 1) is 
uniform, see Example 3.3, so we also have D = 2aU^/^B^^j^^^y^ (n+i)/2' 

Taking n = 2 and s = 1 in (|B.13p we see that the average distance 
between two random points in a circular disc of radius a is Ufa; this is an 
old problem. 

Further examples for a ball of diameter 1 (so a = l/2): ifn=l, then 
ED = 1/3 and ED^ = 1/6; if n = 2, then EL> = 64/457r and ED^ = 1/4; 
if n = 3, then ED = 18/35 and ED^ = 3/10. 

For s = 2, ED^ = 2a^n/{n + 2), as is easily seen directly. 

Example B.IO (Preferential attachment random graph). Pekoz, Rollin 
and Ross [37] have, motivated by the study of vertex degrees in a prefer- 
ential attachment random graph, studied a special case of the triangular 
urn in Section 9 and obtained further results. (They obtain the random 
variable Ka below, for a G {^,1,|,2, ...}, as the limit in distribution, after 
normalization, of the degree of a fixed vertex in one of two slightly different 
random graphs.) 
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In our notation, let Wa, for a > 1/2, be the (limit) variable in Section 
9.1 with a = 2, c = d = 1, Wo = 1 and 6 = 2a — 1. Then (9.1) yields 

ET^^ = r(a) + Res>-1. (B.15) 

r(s/2 + Qj 

(Theorem 2.1 implies that the condition a > 1/2 also is necessary for the 
existence of such a random variable: aG{0,— 1,— 2,...}is clearly impossi- 
ble, since then r(Q) = oo, and otherwise the function in (jB.lSp har p_ < —1 
and /7+ = oo, and if a < 1/2 it is at s = —2a G (/?_, />+), which contradicts 
Theorem 2.1.) 

Pekoz, Rollin and Ross [37] choose a different normalisation, so we define 
Ka ■■= (a/2)i/2w„ and obtain 

In particular, satifies the normalisation EK^ = 1. 

Ka has p+ = oo, p_ = —1 (except when a = 1/2; then /9_ = —2), 
7 = y = 1/2, 5 = 1 - Q, x= iloga, Ci = 2"-V2r(a). 

Pekoz, Rollin and Ross [37] show, among other things, that the random 
variable has the density function 



x„(s) = r{a)J — e-^'/^-u(a -l,^;f-), x > 0, (B.17) 
V avr V 1 la) 



2 . 



where ?7(a, 6; z\ denotes the confluent hypergeometric function of thesecond 
kind; see e.g. [l|, Chapter 13] or [28^] (where it is denoted ^'). This is a con- 
siderably simpler formula than the power series expansion given in Theorem 
9.1. It would be interesting to know whether the density in Theorem 9.1 can 
be expressed using hypergeometric functions also for other triangular urns. 

Note the special case a = 1/2; then ()B.16P simplifies by the duplication 
formula for the Gamma function to 

EK-l/2=T{s/2 + l), Res>-2. (B.18) 

showing that 

= T^'^ = 2-'/^R, (B.19) 

with T ~ Exp(l), see Example 3.2, and R ~ x(2) (the Rayleigh distribu- 
tion), see Examples 3.6 and lB.ll The density function of is thus 

xi/2(x) = 2x6"^", x>0. (B.20) 



Example B.ll (The maximum of i.i.d. exponentials). Let (Ti)'^-^ 
be i.i.d. exponential random variables with Tj ~ Exp(l), and let M„ := 
maxi<j<„rj. Then e~^* ~ U(0, 1), and thus, since e~*^" := mini<j<„ e~^% 

P(e-^^" > x) =P(e-^i > x)" = (1 -x)", 0<x<l, (B.21) 

so e~*^" has the Beta distribution B(l,n). 
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Hence, by Example 3.4, 

Ee^^" = = Res < 1. (B.22) 

Hence M„ has moment generating function of Gamma type. The special 
case n = 1 gives Mi = Ti ~ Exp(l) treated in Example 3.16. 

Mn has /9+ = 1, /9_ = — oo, 7 = 7' = 0, 5 = — n, x = 0, Ci = n!, cf. 
Example 3.4 and Remark 2.8. 

For an alternative proof of ()B.22p . note that if T^^) < • • • < T(„) = M„ 
are Ti, . . . ,T„ arranged in increasing order, then it is a standard observa- 
tion (e.g. by regarding Ti, . . . , r„ as the first points in independent Poisson 
processes) that T(i), T(2) — T(^i), ■ ■ ■ T(^n) £^re independent exponential 

variables with T(^,) — T(fc_i) ~ Exp(l/(n — A; + 1)) (wih T^g) := 0), and hence, 
for Re s < 1 

EgSM„ ^ TT ^ ^ TT i ^ TT _J_ 

LL 1 _ -k + 1) f}^ 1 - si 3 j - s 

= r{n + l) J^]~'^ (B.23) 
r(n + 1 — s) 

Note further that, as n — )• 00, 

EgS(Af„-iogn) ^ r(n + l) _ ^ _ Res < 1, 

r(n+l-s) ^ ^ ^ ^ 



and thus 



(B.24) 

M„-lognAl^, (B.25) 

where W has the Gumbel distribution F{W < x) = which has moment 

generating function Ee*'^ = r(l — s), Res < 1, see Example 3.19. Note 
that it is also easy to prove (1B.25|) directly, since, for a; G M and n large 
enough, 

P(M„ -logn < x) =P(Ti < logn + x)" = (^1 - ^ e"""". (B.26) 

The decomposition 

k=l k=l j=l 

shows also that 



1 

EM„ = J^-, (B.28) 



the n:th harmonic number Hn- We have E = — r'(l) = 7, Euler's gamma, 
and thus, since ()B.24p implies convergence of all moments, 

Hn - log n = E(M„ - log n) ^ E = 7, (B.29) 
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a well-known result by Euler [1^ 

Moreover, it follows from (jR24l) . or from (iB?25]) and (jR29|) . that 

Mn-^Mn^W -¥.W = W (B.30) 

and thus by ([R^ 

oo ^ 

Y.-m-i)^w-^, (B.31) 

where the infinite sum converges in and thus a.s. [l^, Lemma 4.16]. 

Example B.12 (The largest values of i.i.d. exponentials). Generaliz- 
ing Example IB.IH let be the m:th largest of the n i.i.d. exponential 
random variables Ti, . . . ,Tn ~ Exp(l); here 1 < m < n. (The special case 
m = 1 gives M„ treated in Example IB. Ill ) 

Let Ui := e~^' ~ U(0, 1). Then e~*^" ^ is the m:th smallest of the i.i.d. 

uniform Ui, . . . ,Un, and thus e~ " has the Beta distribution B(m,n — 
m + 1). 

Hence, by Example 3.4, 



Thus M^™^ has moment generating function of Gamma type. 

M^™-* has = m, p_ = — oo, 7 = 7' = 0, 5 = — (n — m -|- 1), x = 0, 
Ci = n\/{m — 1)!, cf. Example 3.4 and Remark 2.8. 

Alternatively, (jR32|) can be obtained by the argument in ()B.23p . More- 
over, by the lack of memory for the exponential distribution, — is 
independent of and has the same distribution as M^-i; thus M„ = 

Ml'^^ + Af;^.!, where M;„_ _i is a copy of Mm-i that is independent of Mn"^^; 

this yields Ee"^" = Ee''*^"'"' Ee'^'"-% and follows from ^K22\i . 

As n — 00, 

jEe^CMi-^-iogn) ^ ^-s r(n + l) _ r(m-g) ^ r(m-.) ^ Res < m, 
r(n -|- 1 — s) r(m) r(m) ' ' 

(B.33) 

and thus 

M^'") - log n A W^""^ , (B.34) 

where W^"^^ has the moment generating function of Gamma type 

^ cM/("i) r(m — s) „ N 
Ee'^ — ^, Res<m. B.35 
T[m) 

Comparing with Example 3.1, we see that E e^^'™' = E E"" = E e"" 

and thus W^'^'^ = — logE^, where E^ has a Gamma distribution T{m). 

W^"^^ has /9+ = m, = —00, 7 = 1, 7' = —1, 5 = m — 1/2, x = 0, 
Ci = \p2/K jim — 1)!, cf. Example 3.1 and Remark 2.8. 
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As in (Ib:271) . there is a decomposition 

n— m+l n—m+1 ^ i 

fc=l fc=i i=m 



which shows that 



E ~ ~ ~ Hm-i- (B.37) 



j=m ^ 



Since ()B.33P imphes convergence of all moments, this yields 

jEi^M = E(M^™) - logn) = lim (Hn - H^^i - log n) = 7 - H^^i. 

(B.38) 

Moreover, (IRMI) . (IRM]) and (Ib:38]1 imply 
00 ^ 

^ -(T,- - 1) = - EH^^'") = - 7 + (B.39) 

j=m 



where the infinite sum converges in and thus a.s. 23j, Lemma 4.16]. 

We can also study the joint distribution for several m. In particular, 
(Ib:361) holds jointly for ah m < n, and thus (lR34l) and (lR39ll hold jointly 
for all m. 

Moreover, the conditional distribution of M^™^^'' given Mn \ . . . , 
equals the distribution of the maximum of n — m i.i.d. Exp(l) random vari- 
ables, conditioned on this maximum being at most M,l™\ In particular, 
M^^^ , mP M^"^ form a Markov chain. It is easy to see that this holds 
also in the limit as n — t- 00. Thus W'^^\W^'^\ . . . is a Markov chain, and the 
conditional distribution of T^('"+^) given W'-^\. . . , equals the distri- 

bution of{W\W< W^""^), where W is a Gumbel variabel independent of 
W^'^K Explicitly, for x <y, 

p(^^(m+i) < ^ I ^(m) = = <x\W <y)= exp(e-?^ - e"^). 

(B.40) 

Note that this does not depend on m, so the Markov chain is homogeneous. 
(The chain M^^\ . . . , M^"''' is not.) 



This Markov chain was used by Fristedt [ij] to describe the asymptotic 
distribution of sizes the largest parts in a random partition (after suitable 
normalization); the largest parts have the same asymptotic distribution as 
the largest in a sequence of i.i.d. exponential random variables. 

The Markov chain becomes simpler if we transform to e~^^'"\ The con- 

ditional distribution of e~ " given e~ " , • • • , e~ " equals the distri- 
bution of the minimum e"*"^"-'" oin — m independent uniform random vari- 
ables conditioned on this minimum being at least e~ " ; in the limit it fol- 
lows that the conditional distribution of g-'^''"'''^' given e~^^^\ . . . , g"'^*™' 
equals the distribution of [e~^ \ > e~^'"''), where W' is a copy 
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of W independent of W^'^\ Since e"^' = e'^ ~ Exp(l), it follows that, 
conditionally given 

where T ~ Exp(l) is independent of W^'^\ Consequently, the sequence 
^_vi/(i) ^ ^_vi/(2) ^ ^ ^ ^ ^j^^ same distribution as the sequence of partial sums 
of the i.i.d. Exp(l) sequence Ti, . . . : 

(e-^'^' , e-'^''' ,...) = (Ti , Ti + ,...) . (B.42) 

In particular, this shows again that e"'^''"' ~ T{m) and thus W^'^^ = 
-logTm- 

The same asymptotic distributions Vl^^™-) appear for the largest variables 
in many other situations, see e.g. (27I . Sections 2.2-2.3]. 

Example B.13 (Logistic distribution). Let W have the logistic distri- 
bution with distribution function e^'/(e^' + 1), or equivalently 

— 1 

W>(W > x) = , -00 < X < 00. (B.43) 

+ 1 

By differentiation, the density function is 

1 1 

- = - = - (B 44) 

(e^ + l)2 (e^72 + e-x/2)2 4cosh2(a;/2) ' ^ ' ' 

If Pi has the shifted Pareto distribution with density {x + 1)~^, x > 0, see 
Example 3.14, then P(Pi > e"") = {e-^ + 1)"^ = ¥(W > x) and thus 

I^ = logPi. (B.45) 

As a consequence, by Example 3.14, W has moment generating function of 
Gamma type with 

Ee'^ = EPi" = r(l-s)r(l + s) = -l<Res<l. (B.46) 

smvrs 

Equivalently, W has the characteristic function 

Ee'*^ = r(l-it)r(l + it) = — . (B.47) 

sinh Tit 

We have p+ = 1, /j_ = -1, 7 = 2, 7' = 0, 5 = 1, x = 0, Ci = 27r. 

One way the logistic distribution appears is as the symmetrization of the 
Gumbel distribution. Let W and W be i.i.d. with the Gumbel distribution 
(3.26), see Examples 3.19, and consider W — W, which by (3.35) has the 
moment generating function, for —1 < Res < 1, 

E e^(^-w^') = E e^^ E e''^ = T{1 - s)r(l + s)=Ee'^. (B.48) 
Hence, 

W = W -W'. (B.49) 
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By ()B.49p and (lB.3ip we further have the representation 
1 oo 1 

^ = E 7(^^- - 1) = E -A^^ - ^-^■) 
jyo ■' j=i 

where Tj, j G Z, are i.i.d. with the distribution Exp(l). Since Tj — T-j has 
the moment generating function 

Ee"(^^-^-^) =Ee'^iEe-"^i = — — = -1< Res < 1, 

1 — sl + s 1 — 

(B.51) 

()B.50P is equivalent to 

which is a version of the product formula for sin [l|, 4.3.89] 



oo o 
FT /, 

sm z = z 



3 = 1 ■' 

The random variable with the characteristic function t/ sinht, and thus 
the distribution of VF/vr, is studied by Pitman and Yor [sH] (there denoted 
Si)] among other things, they give the following construction: Let B{t) 
be a standard Brownian motion and let T be the stopping time when an 
independent standard 3-dimensional Brownian motion hits the unit sphere 
in M?. Then 

B{T) = W/-K. (B.54) 

The random variable W /tt (or W , depending on the choice of normaliza- 
tion) appears also as the asymptotic distribution of the rank of a random 
partition, see [§]. 

Example B.14 (Discriminants and Selberg's integral formula). For 

a vector (xi, . . . , x„) of real (or complex) numbers, define 

A(xi, . . . ,a;n) := {xj-Xi). (B.55) 

l<i<j <n 

Thus A(xi, . . . , Xnf' is the discriminant of the monic polynomial with roots 
xi, . . . , Xn- Furthermore, A{xi, . . . , x„) is the well-known value of the Van- 
dermonde determinant '^^^i^^^)" j-i (which apparently was never consid- 
ered by Vandermonde, see [351]). 

Selberg [s^ proved the following integral formula, for n > 2 and Re a > 0, 
Re/3 > 0, Res > max{-l/n, - Re a/(n - l),-Re/3/(n - 1)}, 

/ ••• / \A[xu...,Xn)\^'T\x'^-\l-Xif~Uxi--- dXn 

Jo Jo 
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^ -TT r(a + (j-i).)r(/3 + (j-i).)r(i + j.) 

fii r(a + /3 + (n + j-2)s)r(l + s) " ^ ' ^ 

(For applications of this formula, see e.g. [13] and 0].) This leads to the 
following probabilistic interpretaions, see Lu and Richards [32]. 

For real a,f3 > 0, let Xi,...,Xn be i.i.d. random variables with the 
Beta distribution B(a,/3); then (IB.56P can equivalently be written as the 
expectation 



E\A{Xi,...,Xn[ 



1 2s 



^ -TT r(a + /3)r{a + {j - i)s)r{p + {j - i)s)r{i + js) 

}i r(a)r(/3)r(a + /3 + (n + j-2)s)r(l + s) 

for Res > max{— — Re a/(n — 1), — Re/3/(n — 1)}. This shows that 
A{Xi, . . . ,Xn)'^ has moments of Gamma type. We have 7 = 7' = 0, (5 = 
\ — a — (3 — n/2, = 00 and p_ = max{ — 1/n, —a/{n — 1), —fijiji — 1)} 
(for n > 2). 

Equivalently, ()B.57p shows that the absolute value | A(Xi, . . . , has 
moments of Gamma type. In this case, see Remark 2.8, 7 = 7' = 0, (5 = 
1 — a — /? — n/2, />+ = 00 and p_ = 2max{— —a/{n — 1), —(3/{n — 1)} 
(for n > 2). 

Note that for n = 1, A(Xi) = 1 is trivial, so the simplest non-trivial case 
is n = 2, when (IB.57P says that if Xi,X2 ~ B(a,/3) are independent, then 

^ I ^ ^ i2s _ r(« + fi? r(a + sW + g)r(i + 2s) 
' ^ ^' r(a)r(/3) r(a + /? + s)r(a + /3 + 2s)r(i + s) • ^ • ^ 

We obtain further results by taking suitable limits above, cf. [3|. First, 

note that if -'^j ~ B(a, /3), then l3Xj — ^ Yj ~ r(a) as /3 — )• oo. (For example 
by the method of moments, see (3.6) and (3.1).) By taking limits in ()B.57p . 
using the facts that 

A(/3Xi, . . . , /3X„) = /3"("-i)/2a(Xi, ...,Xn) (B.59) 

and 

'^"" ^^r(/3)"^ ^1 as /3 ^ oo, for every fixed a, (B.60) 
it follows that if Yi, . . . , y„ ~ T(a) are i.i.d., then 

EiA(r.....,r,jp-^n "(°+';-;';>"<;+^'> (b.6i) 

j=2 r(a)r(i + s) 

for Res > max{— 1/n, —a/{n — 1)}. Thus A(yi, . . . ,1^)^ has moments of 
Gamma type, with 7 = 7' = — n, 5 = (n — l)(a — 1/2), p+ = cxd and 
P- = max{ — 1/n, —a/{n — 1)}. In particular, n = 2 yields 

r(Q + s)r(i + 2s) 
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(For a = 1, when 11,12 ~ Exp(l), this is an immediate consequence of 
the fact that \Yi — I2I ~ Exp(l) by the lack of memory in the exponential 
distribution.) 

Secondly, taking /? = a, if Xj ~ B{a,a), then \/^a{Xj — 1/2) — ^ Zj ~ 
N(0, 1) as a 00. By taking limits in (fB757|) . using p3^ - ([B:H0|) and the 
translation invariance 



A{Xi + a, . . . , Xn + a) = A{Xi, . . . ,Xn), (B.63) 
it follows that if Zi, . . . , Z„ ~ N(0, 1) are i.i.d., then 

E |A(Zi, . . . , Z„)|2« = fl Res > -1/n. (B.64) 

t=2 ^(^ + ^) 

(This also follows by letting a — >■ 00 in (IB.61|) . using {Ya — a)/^/a 
N(0, 1), which for integer a is just the central limit theorem for r(l) = 
Exp(l).) Thus A(Zi, . . . , Zn)'^ has moments of Gamma type, with ^ = 7' = 
n{n — l)/2, (5 = 0, /9+ = 00 and p_ = —1. Using the multiplication formula 
(A. 5) for the Gamma function, (IB.64P can be rewritten as 

n n j — 1 

E |A(Zi, . . . , Z„)|2- = (27r)-"("-i)/^(n!)i/2[]i^-^ H 11 ^(^ + Vi) 

j=l j=2 i=l 

i=l l<i<i<n 

The special case n = 2 now just yields 
E \Z, - Z2\^' = ^^^r^ = ^r{s + 1/2), Res > -1/2, (B.66) 



which is immediate because Zi — Z2 ~ N(0, 2), see (3.9). 

The formulas for the moment imply some factorization formulas. Thus, 
a comparison between (|B.64p and Example 3.10 shows the equality in dis- 
tribution 

n 

A(Zi,...,Z„)2 4 (B.67) 

i=2 

where Si/j is stable with index 1/j, and the variables are independent. 
Similarly, ()B.65P and (3.1) show the alternative factorization [33 | 

n 

A{Z,,...,Znf = l[f n (^-68) 

j=l l<i<i<n 

with Gij ~ ^(i/j) independent. 

Similarly, (IKMj) and (3.1) yield 

n 

A(yi, . . . , y„)2 4 A(Zi, . . . , z„)2 J] vf\ (b.69) 

i=2 
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where Vj ~ r(a) are independent of each other and Zi, . . . , by (IB.68P 
this leads to a factorization of A(Y'i, . . . , 1^)^ into independent Gamma vari- 
ables. Another such factorization is given by [s^], where also similar factor- 
izations of A{Xi, . . . , Xn)^ with X ~ B(a, (3) are given for n < 4. 

Example B.15 (Symmetric stable variables). Consider a symmetric 
stable random variable Sa with characteristic function ip{t) := Ee'*'^" = 
e~'*l", where < a < 2. (With this normalization, the Levy measure has 
density c|x|~"~^, where c = (— 2r(— a) cos ^) ^, see e.g. [2l|, Theorem 
3.3].) 

The (locally integrable function) ^, where < Res < 1, has the 
Fourier transform, in distribution sense, Cs|x|~* for a constant Cg given by 



see e.g. 40|, Theorem IV.4.1]; this means that if ip is in the Schwartz class 
S, then 

/OO f'OO 
\t\'-^i^{t)dt = Cs \x\-'^P{x)dx, (B.71) 
-OO J —OO 

where we define the Fourier transform on S by 'tp{x) := e^^^Tp{t) dt. 
(The fact that the Fourier transform is of this type follows by a simple 
homogeneity argument, and the value of Cg then can be found by considering 
the special case a = 2 below.) 

Since ip{t) — )• rapidly as \t\ — )• oo, Sa has a bounded and infinitely differ- 
entiable density function /(x); however, / is not in S. Thus we regularize. 
Let 7/(x) be symmetric and infinitely differentiable with compact support and 
r/(0) = J^^r]{x) dx = 1, and define, for e > 0, "qeix) := e~^7]{x/e), which 
has the Fourier transform %(x) = ff^ex). We then consider the product 
fsix) := f{x)r]s{x) = f{x)r]{ex), whose Fourier transform is / * % = * r/^; 
the function belongs to S, and by applying (|B.7ip with ip = fir and then 
letting e — >■ 0, it follows that (IB.7ip holds with ip{x) = f{x) too, and thus, 
using (iRTOjl 

/CO poo 
\x\-'f{x)dx = c-^ \t\'-^^it)dt 
-OO J —CO 

roo /"OO 

= 2c7^ / f^^e~'" dt = 2c7^a-^ / u^/"-^e-"dn 



2~s^-i/2£i¥)£(i±ll^ (B.72) 

r(i + §) 
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We have proved this for < Res < 1, but by analytic continuation, it 
extends to — a < Re s < 1, and thus 

E|^„r = 2^ Aw-, .n" > -l<Res<Q. (B.73) 

Hence, Sa has moments of Gamma type. We have = a (except when 
a = 2; then /)+ = oo) and /)_ = —1; furthermore, 7 = 1/a, 7' = 1 — 1/a, 
5 = 0, X = log a, Ci = ^J^Ja. 

In the special case a = 2, we have (with our choice of normalization) 

So^ = V2N with N ~ N(0,1), and thus ^7J^ is equivalent to (3.9). (As 
said above, this yields a method to calculate c^.) 

In the special case a = 1, Si has a Cauchy distribution with density 
l/(7r(l + 2;2)), see Example IB. 41 In this case (IB.73|1 yields, using (A. 3) and 
(A.6), 

Vvrr(l-|) vr cos^ 

(B.74) 

This is the same as ()B.6p . and also as (3.13) with n = 1. Indeed, it is 
well-known that Ti = Si, i.e., 71 has a Cauchy distribution. 

Example B.16 (Products of Cauchy variables). Let Xi,X2,... be 
i.i.d. random variables with the Cauchy distribution in Example IB. 41 and 

let Ilfc := Yli X'i be the product of k such variables. (Note that Xi = X~^; 

thus e.g. also 112 = Xi/X2.) It follows from Example IB. 41 that jllfcl has 
moments of Gamma type 

1 

77^^ \2 ' 2J ^ \2 2) cos'=(7rs/2) 

(B.75) 

We have p+ = 1, p_ = -1, 7 = A;, 7' = 0, (5 = 0, x = 0, Ci = 2^. 

The density of Hi = X\ is l/(7r(l + x^)), and the density of 112 = X\X2 

is 

2 log |x| 



Eln.r = + ^1 Vfi - = ,^ -1 < Res < 1. 



7r2(2;2-l)^ 



-cxD < X < 00, (B.76) 



see e.g. Pace [36]. Formulas for the density of 11^ for any integer /c > 1 are 
given by Bourgade, Fujita and Yor [1]. 

Example B.17 (Generalized hyperbolic secant distribution). The 

distribution of the Levy stochastic area A in Example 3.20 is also known 
as the hyperbolic secant distribution, since both the density function and 
the characteristic function are given by the hyperbolic secant 1/cosh (up to 
normalization constants). This distribution is infinitely divisible, and thus, 
there exists a Levy process Ct, t > 0, such that Ci = A; consequently Ct 
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has the characteristic function, cf. (3.36), 



Ee'"'^' = — sGM. (B.77) 
cosh s 



The density is 



7rr(t) 



t + ix 



2 , , ^e^, (B.78) 

see e.g. Pitman and Yor |38l | where many further results are given. 

When t = k is an integer, is the sum of k independent copies of A, so 
by Example 3.20, has moment generating function of Gamma type, with 

E.-&=.-r(l + i)V(l-i)' = ^. |Re»|<f; (B.79) 

\2 IT/ \2 vr/ cos*^ s 2 

we have p± = ±tt/2, 7 = 2k/TT, 7' = 0, 5 = 0, x = 0, Ci = 2^. On the 
other hand, if t is not an integer, then Ct does not have moment generating 
function of Gamma type, since the characteristic function ()B.77p then cannot 
be extended to a meromorphic function in C. 
The density (jRTSl) is for i = 1 



as stated in Example 3.20, and for t 



(B.81) 



2sinhf^ 



Similarly, for every integer t = k > 1, the density (IB.78P is a polynomial in 
X divided by cosh(7rx/2) (k odd) or sinh(7ra;/2) {k even); see Harkness and 
Harkness [la] for explicit formulas. See also [13] for an application. 

Note that = — log ]nfc], where 11^ is the product of Cauchy variables in 
Example lB.161 this is an immediate consequence of the case k = 1 mentioned 
in Example IB. 4i 

As a curiosity, we remark also that the distribution of C2 is related to the 
logistic distribution in Example IB. 131 in the sense that the density function 
of one distribution equals, up to constant factors and a rescaling, the char- 
acteristic function of the other, see (fR44l) . (fR47l) . (lB7f7l> . (IRSTI) . In other 
words, the two density functions are essentially the Fourier transforms of 
each other. 

Example B.18 (Lamperti variables). Let < a < 1 and consider 
La '■= Sa/S'^ where Sa, S'^ are two independent copies of the positive stable 
variable in Example 3.10; thus Ee~*'^° = Ee~*'^" = e"*", t > 0. By (3.16), 
La has moments of Gamma type, using (A. 6), 

r(l-s/a)r(l + s/a) 



T{s/a)r{l — s/a) sin(7rs) 



ar(s)r(l — s) asin(7rs/Q)' 



-a<Res<a. (B.82) 
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We have p± = ita, 7 = 2a^^ — 2, 7' = (5 = x = 0, Ci = 1/a, cf. Remarks 
2.8 and 2.10. 

It is somewhat simpler to consider the power L" = M'^/Ma where Mq, M'^ 
are i.i.d. with the Mittag-Leffler distribution in Example 3.11. By ()B.82p . 
cf. (3.17), 

E(gr = ,;''^"!^i;t'' =^^. -l<Re»<L (B,83) 

T [1 — as)T [1 + as) asm(7rs) 

We now have p± = ±1, 7 = 2 — 2a, 7' = (5 = x = 0, Ci = 1/q, cf. Remark 
2.9. 

The density of L° can be found by Fourier inversion, see e.g. 4l|, p. 445], 
and can be written as 

sin(7ra) 1 

vra + 2 cos(7ra)x + 1 ' 
Consequently, the density of is 

sin(7ra) x°'~^ 



a; > 0. (B.84) 



. . 1 7, 2; > 0. (B.85) 

vr + 2cos(7ra)j;° + 1 ^ ^ 

The random variable was studied (at least implicitl y) by Lamperti 
[2^, and is therefore called a Lamperti variable by James |1^. where also 
further references are given. 

In the special case a = 1/2, (1B.83P simplifies to 1/ cos(7rs/2), so L-^i^ is 
the absolute value of a Cauchy variable, see Example lB.41 which also follows 
directly from ()B.84p . 

Kotz and Ostrovskii (23 | defined, for 0<a</3<2, a random variable 

ya,p = (^a/^)^^" = ^a/^- (The defined F^,/? by giving its density function; 
that the definitions are equivalent follows from (lR84l) .) By (Ib:821) or (iRSSl) . 
Ya^p has moments of Gamma type 

"•^ "//^ r(l + s//3)r(l -s//5) asin(7rs/a)' 

(B.86) 

We have P± = ±a, 7 = 2a~^ — 2, 7' = (5 = x = 0, Ci = 1/q. 

James [18| also considers the more general X^fi ■= Sa/Sa,e, for a > 
and 9 > —a, where Sa and Sa,9 are independent, Sa is a stable variable as 
above, and S^fi has a distribution that is the same stable law tilted by x~^, 
see Remark 2.11. Thus Sa,e has moments of Gamma type given by 

E^r' r(l + g) T{l-s/a + e/a) ^ ^ 

(B.87) 

and X^^g has moments of Gamma type given by, for —a — 9 < Re s < a, 

rc^s^ as r(i + 9) r(i - ./«)r(i + s/a + e/a) 
E X,,, - E 5„ E - ^ ^^^^ r(i - .)r(i + . + e) • ^^-^^^ 
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We have /9+ = a, p_ = -a - 6', 7 = 2(l/a - 1), 7' = 0, (5 = e{l/a - 1), 
X = 0. 

Example B.19 (A generalized exponential distribution). Let /3 > 

and let Vg be a positive random variable with the density function 

TjiTWf"' ' " 

A simple change of variables verifies that this is a probability density func- 
tion, and more generally that, for Res > — 1, 

EV^ - ^ Hx^e-^' dx - ^^^/^ + - ^('/^ + yf^) 

'~ r(i + 1//3) io /3r(i + - r(i//3) ' 

(B.90) 

Vg thus has moments of Gamma type, with p^. = 00, p_ = —1, 7 = 7' = 

Note that /3 = 1 gives the exponential distribution in Example 3.2. In 
general, the distribution of Vp can be seen as a tilted version of the Weibull 
distribution in Example 3.7. 



Example B.20 (Linnik distribution). The Linnik distribution [3l[ | has 
characteristic function 

(B.91) 

where < a < 2. As shown by Devroye 0], a random variable with this 
distribution is easily constructed as 

Xo. := SaVl^'', (B.92) 

where Sa is the symmetric stable random variable in Example IB. 151 ^1 has 
the exponential distribution Exp(l), and these are independent. 
By (iBTfgj) and (3.2), for - min(a, 1) < Res < a, 

E|A^r = E|5^|^Ey//" = 2^ ^ ^ Ml^ + aJ^y a) _ (B^93) 

v^r(i - 1) 

Hence Xa has moments of Gamma type, with p+ = a, p_ = — min(Q,l), 
7 = 2/a, 7' = 1, 5 = 1/2, X = 0, Ci = ^/2^. 

More generally, Devroye [3] showed that if < a < 2 and /3 > 0, and 5q, 
is as above and as in Example IB. 191 and independent of Sa, then 

Xa,fi ■■= SaV^'"" (B.94) 

has characteristic function 

7 rrrrm- (B.95) 

(l + |t|°)V/3 ^ ' 

(This implies that Xa = A^^i, and more generally every Xa^p, is infinitely 

divisible, and that there is a Levy process A^,*, t > 0, such that Xa^t = 
Xa^i/t for all t > 0.) 
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By (Ib:73]) and (IRM]) . for -a/p < Res < a, 



E I X„ « I ^ = E I S„ r E = 2^ ^ 7 ^" . (B.96) 

Hence Xq,.^ has moments of Gamma type, with /)+ = a, /)_ = — a//3, 
7 = 2/a, 7' = 1, 5 = 1//3 - 1/2, X = 0. 

Kotz and Ostrovskii 24] showed that Xa = X^y^^^ where Ya^p is as in 
Example IB . 181 and independent of Xp; this foHows also directly from ()B.86P 
and ()B.93p . For the Linnik distribution see further j25|. 



References 

[1] M. Abramowitz &: I. A. Stegun, eds.. Handbook of Mathematical Func- 
tions. Dover, New York, 1972. 

[2] D. Aldous, Tree-based models for random distribution of mass. J. 
Statist. Phys. 73 (1993), 625-641. 

[3] G. W. Anderson, A. Guionnet, & O. Zeitouni, An Introduction to Ran- 
dom Matrices. Cambridge Univ. Press, Cambridge, 2010. 

[4] P. Bourgade, T. Fujita & M. Yor, Euler's formulae for C{2n) and prod- 
ucts of Cauchy variables. Electron. Comm. Probab. 12 (2007), 73-80. 

[5] M. Bousquet-Melou, Limit laws for embedded trees. Applications to 
the integrated superBrownian excursion. Random Struct. Alg. 29, no. 
4, 475-523 (2006). 

[6] M. Bousquet-Melou &: S. Janson, The density of the ISE and local 
limit laws for embedded trees. Ann. Appl. Probab. 16, no. 3, 1597-1632 
(2006). 

[7] L. Devroye, A note on Linnik's distribution. Statist. Probab. Lett. 9 

(1990), no. 4, 305-306. 
[8] L. Devroye & S. Janson, Distances between pairs of vertices and vertical 

profile in conditioned Galton- Watson trees. Random Struct. Alg. 38 

(2011), no. 4, 381-395. 
[9] P. Diaconis, S. Janson and R. C. Rhoades, Note on a partition limit 

of Bringman-Mahlburg-Rhoades. Bull. London Math. Soc, to appear. 

larXiv:120 5.1252T1 
[10] M. Drmota, Random Trees. SpringerWienNewYork, Vienna, 2009. 
[11] D. Dufresne, The beta product distribution with complex parameters. 

Comm. Statistics - Theory and Methods 39 (2010), no. 5, 837-854. 
[12] L. Euler, De progressionibus harmonicis observationes. Commen- 

tarii academiae scientiarum imperialis Petropolitanae VII (1734-1735; 

printed 1740), 150-161. Reprinted in his Opera Omnia, Series 1, Volume 

14, 87-100. 

[13] P. J. Forrester & S. O. Warnaar, The importance of the Selberg integral. 

Bull. Amer. Math. Soc. 45 (2008), 489-534. 
[14] B. Fristedt, The structure of random partitions of large integers. Trans. 

Amer. Math. Soc. 337 (1993), no. 2, 703-735. 



FURTHER EXAMPLES WITH MOMENTS OF GAMMA TYPE 



19 



[15] J. M. Hammersley, The distribution of distance in a hypersphere. Ann. 
Math. Statist. 21 (1950), no. 3, 447-452. 

[16] W. L. Harkness &: M. L. Harkness, Generalized hyperbolic secant dis- 
tributions. J. Amer. Statist. Assoc. 63 (1968) 329-337. 

[17] L. Hoist, Euler's ({2) = 7r^/6 and the hyperbolic secant distribution. 
Preprint, 2011. 

[18] L. F. James, Lamperti-type laws. Ann. Appl. Probab. 20 (2010), no. 4, 
1303-1340. 

[19] S. Janson, Left and right pathlengths in random binary trees, Algorith- 
mica, 46 (2006), no. 3/4, 419-429. 

[20] S. Janson, Moments of Gamma type and the Brownian supremum pro- 
cess area. Probability Surveys 7 (2010), 1-52. Addendum. Probability 
Surveys 7 (2010), 207-208. 

[21] S. Janson, Stable distributions. Preprint, 2012. arXiv : 1112 . 0220v2 7 

[22] S. Janson, M. Kuba & A. Panholzer, Generalized Stirling permutations, 
families of increasing trees and urn models. J. Combin. Theory Ser. A, 
118 (2010), 94-114. 

[23] O. Kallenberg, Foundations of Modern Probability. 2nd ed.. Springer, 
New York, 2002. 

[24] S. Kotz & I. V. Ostrovskii, A mixture representation of the Linnik 
distribution. Statist. Probab. Lett. 26 (1996), no. 1, 61-64. 

[25] S. Kotz, I. V. Ostrovskii & A. Hayfavi, Analytic and asymptotic prop- 
erties of Linnik's probability densities. I, II. J. Math. Anal. Appl. 193 
(1995), no. 1, 353-371; no. 2, 497-521. 

[26] J. Lamperti, An occupation time theorem for a class of stochastic pro- 
cesses. Trans. Amer. Math. Soc. 88 (1958), 380-387. 

[27] M. R. Leadbetter, G. Lindgren & H. Rootzen, Extremes and Related 
Properties of Random Sequences and Processes. Springer- Verlag, New 
York, 1983. 

[28] N. N. Lebedev, Special Functions and Their Applications. (Translated 
from Russian.) Dover, New York, 1972. 

[29] J.-F. Le Gall, Spatial Branching Processes, Random Snakes and Par- 
tial Differential Equations. Lectures in Mathematics ETH Ziirich, 
Birkhauser, Basel, 1999. 

[30] J.-F. Le Gall & M. Weill, Conditioned Brownian trees. Ann. Inst. H. 
Poincare Probab. Statist. 42 (2006), no. 4, 455-489. 

[31] Yu. V. Linnik, Linear forms and statistical criteria. I, II. (Russian) 
Ukrain. Mat. Zurnal 5, (1953), 207-243; 247-290. English transl. in 
Selected Transl. Math. Statist, and Prob. 3 (1963), 1-40; 41-90. 

[32] I. Lu & D. Richards, Random discriminants. Ann. Statist. 21 (1993), 
no. 4, 1982-2000. 

[33] J.-F. Marckert, The rotation correspondence is asymptotically a dilata- 
tion. Random Struct. Alg. 24 (2004), no. 2, 118-132. 



20 



SVANTE JANSON 



[34] A. M. Mathai, R. K. Saxena & H. J. Haubold. The H-Function. Theory 
and Applications. Springer, New York, 2010. xiv+268 pp. ISBN: 978-1- 
4419-0915-2 

[35] J. J. O'Connor &: E. F. Robertson, Alexandre- Theophile Van- 

dermonde. The MacTutor History of Mathematics archive. 

|http : //www-history . mcs . st-andrews . ac . uk/Biographies/Vandermonde . html 
[36] L. Pace, Probabihstically proving that C(2) = vr^/6. Amer. Math. 

Monthly 118 (2011), no. 7, 641-643. 
[37] E. A. Pekoz, A. Rolhn & N. Ross, Degree asymptotics with 

rates for preferential attachment random graphs. Preprint, 2011. 

'arXiv: 1108.5236. 
[38] J. Pitman & M. Yor, Infinitely divisible laws associated with hyperbolic 

functions. Canad. J. Math. 55 (2003), no. 2, 292-330. 
[39] A. Selberg, Bemerkninger om et multipelt integral. Norsk. Mat. Tidsskr. 

24 (1944), 71-78. 

[40] E. M. Stein and G. Weiss, Introduction to Fourier Analysis on Euclidean 

Spaces. Princeton Univ. Press, Princeton, NJ, 1971. 
[41] V. M. Zolotarev, Mellin-Stieltjes transformations in probability theory. 

Teor. Veroyatnost. i Primenen. 2 (1957), 444-469. 

Department of Mathematics, Uppsala University, PO Box 480, SE-751 06 
Uppsala, Sweden 

E-mail address: svante.janson@math.uu.se 
URL: http : //www. math. uu. se/~svante/ 



